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A gas of self-avoiding surfaces with an arbitrary polynomial coupling to the gaussian curvature and an extrinsic 
curvature term can be realized in a three-dimensional Ising bcc lattice with three local couplings. Similar three 
parameter realizations are valid also in other lattices. The relation between the crumpling transition and the 
roughening is discussed. It turns out that the mean area of these surfaces is proportional to its genus. 



1. Introduction 

It is known that a gas of self-avoiding random 
surfaces[l,2] (SARS) with unconstrained topol- 
ogy and weighted by the usual area term e~^ A 
[3] belongs to the same universality class of the 
Ising model [4, 5]. We discuss the behaviour of this 
gas when other kinds of couplings are introduced , 
like the extrinsic curvature term or a generic poly- 
nomial in the Riemannian curvature. 

There are some three-dimensional lattices 
where the self-avoidance constraint is automat- 
ically satisfied and the SARS are simply the 
boundaries of the spin clusters[6,5,7,8]. In par- 
ticular, in the body centered cubic lattice (bcc), 
it is possible to establish a one-to-one correspon- 
dence between the configurations of Ising spins 
and those of a gas of SARS in the dual lattice 
A in the following way: for each pair of near- 
est neighboring sites < i,j > with opposite spins 
o"j = ±1 , <Ji ^ <Jj we draw the dual plaquette, 
which is a hexagon, and similarly for each pair of 
next to the nearest neighboring sites <C i, J ' of 
opposite spins we draw the dual plaquette which 
now is a square. The special property of the A lat- 
tice is that each link belongs to only three plaque- 
ttes and this allows to enforce the self-avoidance 
constraint. As a consequence, each cluster of like 
spins is bounded by an oriented self-avoiding sur- 
face made with hexagons and squares. For in- 
stance, the surface associated to a cluster of a sin- 
gle site is a polyhedron of 14 faces which is noth- 
ing but the Wigner-Seitz cell of the bcc Bravais 



lattice. Any connected aggregate of these polyhe- 
dra defines a self-avoiding surface. The whole A 
lattice is obtained by filling the space with such 
polyhedra. 

Each potential vertex of a SARS of A is dual to 
a tetrahedron of the bcc lattice selected by a pair 
of orthogonal, adjacent links belonging to the two 
cubic sublattices (see fig. 1) . 
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Figure 1. The tetrahedron 1234 of the bcc lattice 
is dual to a potential vertex of a SARS on the 
dual lattice. The bonds < 1, 2 > and < 3, 4 > 
are dual to squares, while the bonds connecting 
nearest neighboring sites, like for instance 
< 2, 3 > or < 1, 4 >, are dual to hexagons. 



Since these tetrahedra fill the whole space, 
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can write the Hamiltonian of the Ising model as 
the sum of their contributions. On the other 
hand, the spin configurations of each tetrahedron 
can be used to read out the geometrical properties 
of the surface. 

In particular, the Riemannian curvature of 
these polyhedric surfaces is obviously concen- 
trated on the vertices. The contribution R v to 
the curvature of the vertex v dual to the tetrahe- 
dron drawn in fig.l can be written as a function 
of the spin configuration of the four sites 1,2,3,4: 

24 R v = <7i3 + <7i4 + <7 23 + <T24 - 3<7i 2 34 ~ 1 (1) 

where <7 8 j = <Ji<Jj. Note that when the four 
vertices 1,2,3,4 belong to the same cluster then 
R v = as it should be. 

Denoting by S a generic, not necessarily con- 
nected SARS, by N S (S) and Nh(S) the number 
of squares and hexagons belonging to S, and by 
x(S) its Euler characteristic, where 

x(S) = disconnected comp. — ^handles) (2) 

the partition function Z(j3 s , f3h, fit) of a gas of 
SARS can be written as 



j-i3 s N s (S)-l3' h N h (S)-l3 tX (S)] 
S 



(3) 



where, besides the area terms [3 S and fth there is 
also a "topological fugacity" e~^ 1 . 
By using the Gauss-Bonnet theorem 



J d 2 i^R=Y,Rv=x(S) 



(4) 



and the previous expression for R v , such a parti- 
tion function can be rewritten as an Ising model 
on the bcc lattice: 



(7 = ±1 



(5) 



with the (reduced) Hamiltonian given by 

<ij> v 

where a v = <7 8 '<7j<7/;<7/ is the product of the spins 
of the four vertices of the tetrahedron dual to the 
potential vertex v, fih = fi'h ~ ano - an uninter- 
esting additive constant has been omitted. 



2. Riemannian Curvature 

A crucial point is that any power of the curva- 
ture R v can be expressed again as a suitable com- 
bination of the three kinds of couplings < ij >, 
<C hi ^> and (ijkl), as it follows directly from 
eq. (1). Thus any partition function of SARS, 
weighted by an action A polynomial in the cur- 
vature R 



A 



2 C^/g[fi + /3 t R + aR 2 



(7) 



where /j, is the cosmologic constant, can be dis- 
cretized as an Ising model on the bcc lattice de- 
scribed by a partition function Z(f3 s , fih, fit) of 
the same functional form of eq.(5), where ft s ,f3h 
and /3t are suitable functions of the parameters 
H, fit, a, . . .. For instance, a term proportional to 
R 2 has been introduced in some models [9] of dy- 
namically triangulated surfaces (DTS) in order to 
study the transition from the crumpled to the col- 
lapsed phase, which is favoured for a < 0. If we 
add a R 2 term to eq.(6) the new couplings are 



13 s = j3 s 



Sit = Sit 
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(8) 



(9) 



(10) 



The constraint of self-avoidance forbids collaps- 
ing of our surfaces. Note however that there are 
values of a < and fit such that (i^ = and 
fit = where the lattice decouples into a pair 
of simple cubic lattices. In these sublattices self- 
intersecting, collapsed surfaces are allowed. Note 
that the numerical coefficients of the above equa- 
tions do not have any universal meaning, being 
of course specific of the bcc lattice. It is however 
worth noting that also in the other lattices where 
the self-avoidance constraint is automatically sat- 
isfied, one can describe a gas SARS coupled to an 
arbitrary polynomial of the curvature R in term 
of an Ising model with only three kinds of cou- 
plings. This is obvious in the lattice defined in 
[6], once one realizes that it is topologically dual 
to a bcc lattice in a disguised form. It is also 
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true in the dual lattice C of the face centered cu- 
bic (fee) lattice, used in ref.[5], where apparently 
there are more complicated couplings. Indeed, 
there are two kinds of potential vertices for SARS 
of £ : one is specified by the spins of four suitable 
neighboring sites, while the other needs the spin 
status of eight sites. Denoting respectively by i?' 
and R° the contribution to the curvature of these 
vertices and using formulas similar to eq.(l), one 
easily verifies that 



and 



(K? = -^K 



(ID 



(12) 



Hence, a gas of SARS described by the action (7) 
may also be discretized as a Ising model in the fee 
lattice with three kinds of couplings associated to 
R\j, R° and to the pairs of nearest neighboring 
sites. 

3. Extrinsic Curvature 

Coming back to the bee lattice, we want to show 
that it is also possible to put in the same three- 
parameter form a SARS partition function with 
an extrinsic curvature term 



E = ^(i-C 08 ^^) 

(a,b) 



(13) 



where {} a b denotes the angle between two incident 
faces a and b . In the A lattice there are only two 
kinds of such angles: if a and b are both hexagons 
one finds cos(^ a 5) = 



square, then cos(^ a 5) 



i ; if one of them is a 



Modifying the partition function Z(f3 s , fih, fit) 
of eq.(6) as follows 



Z E {l3,J h J t ) = J2 



-/3H — XE 



(14) 



<7 = ±1 



where e A is the fugacity associated to the ex- 
trinsic curvature term, we find 

A 



Ph = Ph + 



(15) 



p s = (] s + cX 



32^-21 



24 



Pt = Pt 



11A 
~6~ 



(16) 



(17) 



In the DTS models the coupling to E is used 
to control the crumpling transition [10] : A drives 
the system from a crumpled phase to a smooth 
phase. We see from the above expressions that 
the flow towards the smooth surfaces implies a 
cooling of the gas. It is conceivable that the crum- 
pling transition of the single surface would corre- 
spond to the roughening point of the Ising model. 
A numerical analysis in this direction might be 
very interesting. 

4. Condensation of Handles 

Microscopic analysis of the topology of the 
SARS near the critical point of the Ising model 
by directly counting the number of handles of 
the connected boundaries of the spin clusters has 
been performed only recently [8,11]. 

It has been found for different lattices and in a 
wide range of temperatures near the critical point 
of the Ising model, that the number Af of SARS 
as a function of the area A and the genus G is 
accurately described for large surfaces by [8,11] 

Af(xA f ^e-^ A (18) 

where the "cosmological constant" /j, is a function 
of the lattice shape and f(G) is an almost linear 
function of the genus. It has been pointed out 
[11] that SARS with unconstrained topology can 
also be reconstructed from the spin configurations 
of a simple cubic lattice, even if it is well known 
that in such a lattice the cluster boundaries do 
not satisfy the self-avoiding constraint, because 
of the possible presence of lines of self-contact. 
These lines are the only potential sources of am- 
biguities in the topological reconstruction of the 
surface. By small deformations of these surfaces 
it is possible to eliminate, at least locally, these 
contact lines. However, if these deformations are 
chosen arbitrarily it may appear global obstruc- 
tions in the surface reconstruction. It has been 
found a set of rules to disentangle in a simple way 
these contact lines which avoids global obstruc- 
tions and assigns to each Ising configuration a set 
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of SARS [11]. Notice that this set of rules is by no 
means unique: there are other sets of consistent 
rules yielding for the same Ising configuration a 
different set of SARS. Notwithstanding, the only 
modification observed in eq.(18) is a variation of 
the parameter ji. 
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Figure 2. The mean area of an interface as a 
function of its genus in a cubic lattice. See ref. 
[11]. 

An important consequence of the linearity of 
f(G) is the phenomenon of condensation of han- 
dles: the mean area < A >a of the surface at 
fixed genus is proportional to its genus. Indeed, 
by eq.(18) 



< A > ( 



9 i K r 



one gets at once 



< A > ( 



f(G) 



-G ■ 



(19) 



(20) 



[i [i [i 

Thus the number of handles per unit area is a 
non-vanishing constant b//i; this means that the 
SARS are dominated by configurations with a 
large number of microscopic handles, then the 
cluster boundary is much more similar to a sponge 
rather than to a smooth surface. See fig. 2. 

The universal parameter b is related to the 
string susceptibility j s by b = 2 — j s . The fact 
that all SARS in the bcc lattice can be generated 
by gluing together in all the possible ways the 



Wigner-Seitz cells suggests that they belong to 
the same universality class of the branched poly- 
mers, which gives f s = ^. Preliminary numerical 
simulations suggest a slightly larger value. 

There are arguments suggesting that SARS 
with unrestricted topology could be in a different 
universality class [12,13]; however recent numeri- 
cal simulations on self- avoiding, flexible, dynam- 
ically triangulated surfaces seem to indicate that 
the number of handles is irrelevant when classi- 
fying SARS in universality classes [14]. It would 
be very interesting to evaluate the radius of gy- 
ration of SARS in the dual of bcc lattice in order 
to settle this point. 
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